ohi logo
OHI-Northeast | OHI Science | Citation policy

1 Summary

3 Data Cleaning

## # A tibble: 6 x 6
##    year stat_area species                      pounds stock_id stock       
##   <dbl>     <dbl> <chr>                         <dbl> <chr>    <chr>       
## 1  1996         0 CONFIDENTIAL SPECIES       11517161 <NA>     <NA>        
## 2  1996       462 COD                           11827 CODGMSS  GOM Cod     
## 3  1996       462 CUSK                           1065 <NA>     <NA>        
## 4  1996       462 FLOUNDER, AMERICAN PLAICE…    21191 PLAGMMA  Plaice      
## 5  1996       462 FLOUNDER, WITCH / GRAY SO…     5496 WITGMMA  Witch Floun…
## 6  1996       462 HADDOCK                         296 HADGM    GOM Haddock

5 Catch per region

Let’s filter the catch data to just the statistical areas in our region.

## # A tibble: 20 x 12
##     year stat_area species pounds stock_id stock rgn_name rgn_id state
##    <dbl>     <dbl> <chr>    <dbl> <chr>    <chr> <fct>     <int> <fct>
##  1  1996       464 COD      31842 CODGMSS  GOM … Gulf of…      3 <NA> 
##  2  1996       464 CRAB, … 200310 <NA>     <NA>  Gulf of…      3 <NA> 
##  3  1996       464 CUSK     44111 <NA>     <NA>  Gulf of…      3 <NA> 
##  4  1996       464 FLOUND…  13415 PLAGMMA  Plai… Gulf of…      3 <NA> 
##  5  1996       464 FLOUND…  22237 WITGMMA  Witc… Gulf of…      3 <NA> 
##  6  1996       464 FLOUND…   2930 YELCCGM  CC/G… Gulf of…      3 <NA> 
##  7  1996       464 HADDOCK  12570 HADGM    GOM … Gulf of…      3 <NA> 
##  8  1996       464 HAKE, …  55368 HKSGMNGB Nort… Gulf of…      3 <NA> 
##  9  1996       464 HAKE, …   4256 <NA>     <NA>  Gulf of…      3 <NA> 
## 10  1996       464 HAKE, …  25578 HKWGMMA  Whit… Gulf of…      3 <NA> 
## 11  1996       464 LOBSTE… 439729 <NA>     <NA>  Gulf of…      3 <NA> 
## 12  1996       464 MONKFI…  84788 <NA>     <NA>  Gulf of…      3 <NA> 
## 13  1996       464 POLLOCK  30408 POKGMASS Poll… Gulf of…      3 <NA> 
## 14  1996       464 REDFIS…   3960 REDGMGB… Redf… Gulf of…      3 <NA> 
## 15  1996       464 SCALLO…    427 <NA>     <NA>  Gulf of…      3 <NA> 
## 16  1996       464 SKATE,…   2104 <NA>     <NA>  Gulf of…      3 <NA> 
## 17  1996       464 WOLFFI…   1481 WOLGMMA  Wolf… Gulf of…      3 <NA> 
## 18  1996       465 COD      41754 CODGMSS  GOM … Gulf of…      3 <NA> 
## 19  1996       465 CRAB, … 107500 <NA>     <NA>  Gulf of…      3 <NA> 
## 20  1996       465 CUSK     10134 <NA>     <NA>  Gulf of…      3 <NA> 
## # … with 3 more variables: ohi_rgn_prop_area <dbl>, geometry <POLYGON
## #   [m]>, catch <dbl>

Map one species

Visualize the data

## # A tibble: 30 x 7
## # Groups:   species, stock, stock_id, year, rgn_name [30]
##    species stock stock_id  year rgn_name                    rgn_id    catch
##    <chr>   <chr> <chr>    <dbl> <fct>                        <int>    <dbl>
##  1 ALEWIFE <NA>  <NA>      1998 Gulf of Maine                    3 2627.   
##  2 ALEWIFE <NA>  <NA>      1998 Maine                            6  539.   
##  3 ALEWIFE <NA>  <NA>      1998 Massachusetts-Gulf of Maine      7    7.07 
##  4 ALEWIFE <NA>  <NA>      1998 New Hampshire                    9   47.8  
##  5 ALEWIFE <NA>  <NA>      1999 Gulf of Maine                    3   43.4  
##  6 ALEWIFE <NA>  <NA>      1999 Maine                            6    8.90 
##  7 ALEWIFE <NA>  <NA>      1999 Massachusetts-Gulf of Maine      7    0.117
##  8 ALEWIFE <NA>  <NA>      1999 New Hampshire                    9    0.790
##  9 ALEWIFE <NA>  <NA>      2000 Gulf of Maine                    3  241.   
## 10 ALEWIFE <NA>  <NA>      2000 Maine                            6   49.4  
## # … with 20 more rows

Create maps of catch for all species using just the most recent year

Calculate species catch per OHI region

Let’s look at total regional catch for each species (not stock)

Clearly atlantic herring is making up the majority of catch!


There are 155. Way more than we have stock information for. Let’s see what species also have stock ids

## # A tibble: 28 x 3
##    species                                 stock                  stock_id
##    <chr>                                   <chr>                  <chr>   
##  1 COD                                     GB Cod East            CODGBE  
##  2 COD                                     GB Cod West            CODGBW  
##  3 COD                                     GOM Cod                CODGMSS 
##  4 CONFIDENTIAL SPECIES                    Southern Whiting       HKSOGMMA
##  5 FLOUNDER, AMERICAN PLAICE /DAB          Plaice                 PLAGMMA 
##  6 FLOUNDER, SAND-DAB / WINDOWPANE / BRILL Southern Windowpane    FLDSNEMA
##  7 FLOUNDER, SAND-DAB / WINDOWPANE / BRILL Northern Windowpane    FLGMGBSS
##  8 FLOUNDER, WINTER / BLACKBACK            GB Winter Flounder     FLWGB   
##  9 FLOUNDER, WINTER / BLACKBACK            GOM Winter Flounder    FLWGMSS 
## 10 FLOUNDER, WINTER / BLACKBACK            SNE/MA Winter Flounder FLWSNEMA
## # … with 18 more rows


We have records for 0 catch. But what about the missing data? Let’s look at ALEWIFE.

##  [1] 1998 1999 2000 2003 2006 2007 2012 2013 2014 2015 2016 2017

Ok clearly we are missing 2001, 2002, 04-05, 2008-11.

6 Gapfill

We need to gapfill this missing data. When a species/state combination has missing data for a year, we can not assume it has a catch of 0. Since we calculate a rolling average of catch, NAs will remain as NA’s and the average will rely on just one or two years of catch. This is done to account for any wild fluctuations in catch year to year.

7 Remove non-food catch

Some of these species are harvested for food as well as other markets like pet food or bait. We want to make sure this goal captures catch meant for human consumption. We have data from NOAA that identifies the amount of catch per species, state and year meant for food, bait, and other markets. This data was cleaned in prop_catch_food_bait.Rmd.